Parametrization of the QCD coupling in the Evolution Equations 

B.I. Ermolaev 

H. Niewodniczanski Nuclear Physics Institute PAN, 31-342 Krakow, Poland 
and loffe Physico-Technical Institute, 194021 St. Petersburg, Russia* 

S.I. Troyan 

St. Petersburg Institute of Nuclear Physics, 188300 Gatchina, Russia 

We examine the parametrization of the QCD coupling in the Evolution Equations, including 
QQ ' DGLAP. Our conclusion is that the well-known parametrization, where the argument of the coupling 

, is k\/l3 or just k\, stands only if the lowest integration limit in the transverse momentum space 

■ (the starting point /i^ of the -evolution) obeys the relation ^ ^ ^qcd exp (7r/2), otherwise the 

' coupling should be replaced by the more complicated expression presented in Eq. p7|l . 

PACS numbers: 12.38.Cy 
• I. INTRODUCTION 

^ ! 

' By various reasons, parametrization of as in QCD processes has always been among the QCD topics on discussion. 
' (— I 'I In particular, the early summaries on the situation in the hard processes were presented in Refs. (T], Q and one of 
the latest summaries is Ref. [31. It is well-known (see e.g. Refs. [l|,0) that in the Feynman graphs contributing to 
' scattering amplitudes in the hard kinematics, as well as to the Sudakov form-factors, the parametrization of as is 
universal and simple: 
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as{kl), (1) 



I \ with k± being the transverse momentum of the soft gluon momenta. This parametrization is also used in the DGLAP 
J> ' equations in the integral form and leads to the parametrization as = as{Q'^) in the differential DGLAP equations. 
OO Besides this example, there are other cases where it is important to know the parametrization of a^. For example, 
. the parton distributions F and some other objects are often obtained through composing and solving the evolution 
■ equations of the Bethe-Salpeter (BS) type, with one intermediate gluon factorized. Naturally, such equations can be 
, different and can be derived by different ways. In particular, both DGLAP and BFKL equations belong to this group. 
I/"} We do not consider obtaining those equations in the present paper. Instead, our goal is to consider the treatment of 
as in BS equations in general. The parametrization of Us in the BS equations was discussed in Ref. [J]. In the much 
more detailed way the treatment of in the BS equations was suggested in Ref. ^] where the parametrization 



as^as{-kl/P), (2) 

k> , with (3 being the longitudinal moment fraction (see Eq. ([5])), was derived. In principle, the parametrization ^ looks 
^1 natural and well-expected from the preceeding estimates, however the method of Ref. [3| includes serious contradictions 
^ which make the derivation of Eq. ^ unreliable. In the present paper we consider the parametrization of as in the 
Bethe-Salpeter equations. We show what was the mistake done in Ref. [3] and correct it. As a result, we obtain 
a more complicated expression for the effective coupling in the BS equations. Then we demonstrate that when 
the lowest limit of integration over k\ is large enough, the effective coupling can be simplified down to as{k'j_/ (3). 
Earlier, in Ref. we argued that one of the easiest ways of treating as in the BS equation could be the use of the 
Mellin transform. Nevertheless by various reasons, for example for Monte-Carlo simulations^ and for performing 
the numerical calculations in the way convenient for experimentalists^, one may prefer not to involve this operation. 
To this end, in the present paper we examine the treatment of as in the Bethe-Salpeter equations without using the 
Mellin transform. The paper is organized as follows: in Sect. II we remind how the parametrization of as in the 
hard kinematics, i.e. Eq. can be obtained. We suggest in Sect. Ill the way to derive the parametrization of as 
in the BS equations for the scattering amplitudes in the forward kinematics. Due to the Optical theorem the parton 
distributions are proportional to the imaginary parts of the forward amplitudes. It allows us to fix in Sect. IV the as 
-parametrization in the BS equations for the parton distributions. Finally, Sec. V is for concluding remarks. 



* Permanent address 



2 



II. PARAMETRIZATION OF a, FOR QCD PROCESSES IN THE HARD KINEMATICS 

From the technical point of view, calculation of QCD processes in the hard kinematics is much easier compared to 
the Regge kinematics, so to begin with we remind known results for the parametrization of for the hard processes. 
Let us consider the contribution of the Feynman graph depicted in Fig. 1. It correspond to a certain factorization 
of the soft virtual gluon with momentum k propagating in the t -channel. The cases with u and s -channel gluons 
factorized can be considered quite similarly. In the present paper we do not discuss the factorization but focus on 
the treatment of as only. The solid lines in Fig. 1 denote quarks, though the generalization to the case of gluons is 
obvious. Through the paper we will assume that the lower particles, with momenta pi, p'l, have small virtualities ^ /i^ 
whereas virtualities of the upper partons, with momenta q, q' are large; ~q^ ^ ~q'^ ~ ^ /i^. This assumption 
allows us to extend our analysis, with minor changes, to DIS and DVCS. Applying the Feynman rules to Fig. 1 we 
obtain: 

Factor 2s in Eq. ^ appears from simplifying the numerator. We have denoted s = (q+p)"^ and assume that ^/s :§> ^ 
and any of the involved masses. We dropped the color factors in Eq. ^ as unessential for our analysis. M corresponds 
to the blob in Fig. 1. In terms of the Sudakov variables 

k = -a{q + xp) + (3p + kj_ (4) 

and therefore 

k'^ = —safi — k\^ 2pk — —sa, 2qk ~ s{l3 — xa). (5) 

In Eq. ([5]) we have neglected the virtuality = of the initial parton and denoted s = 2pq, x — Q'^ /2pq, = —q^- 
In terms of the Sudakov variables Eq. ([3]) is 



* *4^2 j aaapa^x _ q2 _ _ ^2 + [^^ _ _ (^^^ + ^2 _ • 

Amplitude M is unknown, so it is impossible to perform the integration over any of the variables in Eq. ([6]). However, 
if we assume the leading logarithmic (LL) accuracy, we can use the QCD -generalization [8] of the bremsstrahlung 
Gribov theoremQ. According to it, M does not depend on a and (3. Integrating over a in Eq. © is conventionally 
performed with closing the integration contour down and taking the residue at sa — (+fc^ — ~ It converts 

Eq. IH) into 



1 f\j^-^)..f ^2 ,2^ f k\ 



where Po — x + k'^/s. Obviously, Po ^ x — /2pq when x ~ 1 and the upper limit s of the integration over k'j_ can 
be changed for Q^. The minus sign of the as -argument in Eq. ([7]) indicates explicitly that the argument is space-like 
and for the space-like argument as is given by the well-known expression: 

ai-TT^)- J (8) 

^ (1-/3)^ &ln(fc^/((l-/?)A2)j 

We have used here the standard notation b — 127r/[llA^ — 2n/] and denoted A = Aqcd- Usually the argument 
of as in the Ihs of Eq. ([8]) is written without the minus sign as as(<?^) = l/[b\n{q'^ / A'^)] but this the expression is 
complemented by the statement that q^ < 0. In our case — /3) > 0, so we keep the minus sign to show explicitly 

that the argument of as is space-like. Obviously, Qas = in Eq. ([S]). In contrast to it, when the argument of as is 
time-like, 7^ and in this case 

,9, 1 1 Infm^/A^) + tTT 

' 6[ln(m2/A2) -ztt] 6 [In(m2/A2) -f 7r2] ' ^ ' 

When the contributions ~ tt in Eq. Q are neglected, there is no difference between as(m2) and as(— m2). With 
the LL accuracy, k'j^/(l — (3) « k'^ and therefore in Eq. ([7]) as ~ as(fc^). The minus sign of the argument of as is 
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traditionally dropped, which drives us back to the standard expression Eq. ([T]). Finally, differentiation over In(Q^) 
converts Eq. ([7]) into the well-known form: 

dHQyJ^)^'^L ^^^"^^^'^ ^"4(1-5)) "-^^1 ^M(.,Q,fc,). (10) 

On the contrary, with the leading logarithmic accuracy, (3q k'^ / s a.t x <^ 1 and the integration over fc^ in Eq. ([7|) 
runs from /z^ to s. 

III. PARAMETRIZATION OF a, IN THE BETHE-SALPETER EQUATIONS 

In this section we study the parametrization of in the Bethe- Salpeter equation for the forward scattering 
amplitude A . Let us assume that A obeys the following Bethe-Salpeter equation: 



A = Ao + —^ / dadl3dkiM{{q + ky,Q',{sap + ki))- ^ , ' ' ^, ^\ (11) 



sfc^ as{sa{l - (3) - k]_) 



Ao + -^ [ dkidpdm^M{sp,QMm^f3 + kl))- ~ ^^^^ "^(""'^ 



The second term in the rhs of Eq. (fTTj) is depicted in Fig. 2. Following Ref. [3], we have replaced the Sudakov variable 
a by the new variable m? = [p — k)"^: 

^2 I 1,2 

(12) 



1-/3 



M in Eq. (fTTj) denotes the upper blob in Fig. 2. It includes both the off-shell amplitude A and a kernel. We will 
define M later. Now we just notice that Eq. pT|) can be solved only after M has been known, stands for an 
inhomogeneous term. We do not specify Aq because it does not affect the ois -parametrization. In the first place we 
focus on integrating over a in Eq. and introduce 



/= / dm'M{sl3,Q\ {m'l3 + ki)) ^ ' ^ ,7 ^ „ ^ ^ ^ (13) 



(l-/3)fci a,(m^) 
{•ni?j3 + fcj^ — ie)^ [m? + ie] 

The integrand of Eq. (jl3[) has the singularities in rn^. First, there are two poles from the propagators: 

■ni^ ^ -k\/ P + le (14) 

and 

m^^O-ie. (15) 

Second, there are two cuts. The first cut is originated by the fc^ -dependence of M. In particular, it can be the 
logarithmic dependence. The cut begins at 

= -kl/(3 + ie (16) 

and goes to the left. The second cut is related to a^. It begins at 

m^^O-ie (17) 

and goes to the right. The singularities (|14m7p are depicted in Fig. 3. The integration over in Eq. ([TH]) runs 
along the 5Rm^ -axis from — oo to oo, so the integral can be calculated with choosing an appropriate closed integration 
contour C and taking residues. The contour C should include the line — oo < < oo and a semi-circle with 
radius R. The contour Cu may be situated either in the upper or in the right semi-plane of the -plane. However, 
if we choose Cr to be in the upper semi-plane, we should deal with the cut (|16p of an unknown amplitude M, which 
is impossible without making assumptions about M. Such a closing of the contour was chosen in Ref. [1] but the 
contribution of the cut (|16p was not taken into account because there was made the assumption that 

M{sl3,Q^,e)^ M{sP,Q\kl), (18) 
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i.e. that k'j_ ^ |m^/3|. This assumption made possible to calculate the integral / with taking the residue at the pole 
Eq. (|16p where fc^ = |m^/3|. This obvious contradiction between Eqs. and makes the method of Ref. Q 
inconsistent. Alternatively, choosing the contour Cji in the lower semi-plane involves analysis of the cut (|17p of 
and as is known. By this reason, we choose the latter option for Cr. So, as shown in Fig. 4, the closed contour C 
includes the line — c» < rri^ < oo, the lower semi-circle Cr and the contour Cent which runs along both sides of the 
cut (fT7|) . According to the Cauchy theorem, 

, .2^^ (l-/3)fci asim') 



-27ri-^i-2^M(s/3, -^1/(1 - |3))as{^l^). 



"J. 



The rhs of Eq. (|19p is the residue at the pole and is introduced to regulate the IR singularity for as . It should 
be chosen as large as 

M >> A (20) 

to guarantee applicability of the perturbative expression in Eq. ^ for as- When the initial partons in Eq. ([50)1 are 
quarks, ^ should also obey /i ^ the quark mass. Obviously, 

IC=I + lent + Ir (21) 

where / is defined in Eq. (|13p , Ir stands for the integration over the lower semi-circle and Icut refers to the integration 
along the cut (fTT]) . Ir can be dropped because /_r ^ when oo. Now we specify Icut- 

Icut = -2^ rdmHl{s(3,Q^ (m' p + kl)) ^^WtS^^^^^ ' ^^2) 
J ^2 [m'^p -I- K]_ — ley 

The integration in Eq. (|22p cannot be done precisely because it involves the unknown amplitude M depending on 
771^. Contrary to the case of the hard kinematics considered in Sect. Ill, the arguments of the Gribov bremsstrahlung 
theorem 9] cannot be used here. Nevertheless, it is possible to estimate Icut- Indeed, the m^- dependence of M in 
Eq. ((22|) can be neglected in the region 



< kl/(3- (23) 

Doing so, we obtain the following estimate of Icut'- 

kj_ J ^2 nfi b kj^ ^ J ^2 m? \[iy {m? / K^) + tt^] 

(24) 

When n is chosen as large that 

ln(^/A) > 7r/2, (25) 
we can drop ir"^ in Eq. (|24p and arrive at the estimate 

/c„t« ^^%^M(s/3,Q',fci)[a.(fci//3)-as(M')]- (26) 

Generally, when ^ is not kept very large, though the condition /i ^ A must be satisfied, the integration in Eq. (j24p 
issues that ^ 



^ A similar estimate was done in Ref. ([5|). 
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Obviously, Eq. converts into when fi obeys ([25]). Combining Eqs. (119121^ 1 and ^7^, we rewrite Eq. dTS]) as 
follows: 



^M(.A Q^ - P) [a.i,^) - - arctan (^^"^^l^^^ • (28) 



Assuming that obeys Eq. ((25)) allows to simplify Eq. ((28)) down to the very simple expression free of the infrared- 
dependent contributions: 

— 2?7r 

/ ~ -^M(s/3, Q',kl)a,{kl/P). (29) 



IV. PARAMETRIZATION OF IN EQUATIONS FOR THE PARTON DISTRIBUTIONS 

Now we can easily consider to the Bethe-Salpeter equations for the parton distributions. Indeed, the parton 
distributions are proportional to So, combining Eqs. ((TT)) and ((28)) and taking ^A, we arrive at the following 
expression; 



=^llt i " + h(^^) - ^ (^^ + ln(.i/M^) ln(,7A^) 



^[ln(fci//3A2)-ln(AiVA2)] 



(30) 



In Eq. ((50)) we have used fi^ as the infrared cut-off and the starting point of the integration over k'j_ . 3M 7^ when 
{q + k)^ > 0. In terms of the Sudakov variables the condition s > [q + k)"^ > can be approximately rewritten as 

1> P> Po^x + kl/s. (31) 

It defines the limits of the integration over /3. It also fixes the upper limit of integration over fc^ in Eq. ((50)) . We 
consider below the large- a; and small- a; situations. 



A. Parametrization of as in DGLAP 

When a; < 1, it follows from Eq. ((5T)) that /3o ~ x, so we can write Eq. ((50)) in the DGLAP-like form: 



. .7r2 + ln(/s2^//3A2)ln(/iVA2). 

and express the off-shell 5Af through the off-shell ^A: 

3M((9 + k f, Q\ k\) - 3A((g -f fc)2, fci)i. (33) 

Assuming that obeys Eq. ((^5)) . allows to approximate the expression in the square brackets in Eq. ((5^ by the 
expression in Eq. ((M)) . After that we obtain: 



/ ^(l-m^((9 + fc)',Q^fcl)«.(fci//3)- (34) 

Approximating as(k\l(i') ~ Q;s(fc^) and adding non-ladder contributions, where the argument of is given by 
Eq. ([7]), leads to multiplying the rhs of Eq. ((55)) by the LO DGLAP splitting function(s) which we denote P(/3) 
without specifying. After that, differentiating with respect to j [i}) converts Eq. ((5^ into the well-known 

DGLAP equation: 

f'<'^.p,^,a^„^,<3,,. ,35, 



91n(Q2/^2) 27r }^ (3 

Let us remind that, in general, the parametrization of for the contributions of the non-ladder graphs (see Eq. ((7))) 
considerably differs from the parametrization of the ladder contributions in Eqs. ()30l34p . 
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B. Parametrization of 0:3 in small-a; evolution equations 

When a; <C 1, the lowest Umit of integration in Eq. (|30p is s/3o « fc^. Also the upper limit for the /? and fc^ 
-integrations is s/3 ~ 1. Besides, there can be a kernel K which should not be associated with the DGLAP splitting 
functions. So, in the small- a; limit Eq. ([50)1 can be written as follows: 

However, the change of the integration limits does not affect the parametrization of a.,- Indeed, if /i obeys Eq. (j25p . 
the expression in the squared brackets in Eq. (j36p can again be simplified down to as{k'^/ P), otherwise it remains as it 
is in Eq. (|36p . In contrast to Eq. , the small- a; Eq. ([55)1 cannot be simplified down to Eq. ([55)1 with differentiating. 

V. DISCUSSION 

The analysis of the parametrization of as we have done in the present paper can be addressed the wide group of 
existing and forthcoming evolution equations of the Bethe-Salpeter type, including BFKL and DGLAP, where one 
virtual gluon is factorized out of the blob. Such a gluon can propagate in the s-channel as well as in the crossing 
channels. We demonstrated that basically the parametrization of as depends on the channel. For the crossing channels, 
where the factorized gluons are soft, the parametrization of as is universally given by Eq. ([Tj): as = as{k'^/{l — /?)). 
The case of the s -channel, where the factorized gluon is not soft, is more involved. The effective coupling ag^^ here 
is given by the following expression: 

^ ' nb V7r2 -f ln(fc2^//3A2)ln(^2/jY2)y ^ ^ 

However, if the lowest limit /z^ of integration over k'j^ is chosen large enough to obey 

» Aqcdc^/^ ~ 5Aqcd, (38) 

i.e. when /i > 50Aqcd, Eq- (EH) can be simphfied: 

aff « as{kl/l3). (39) 

We obtained Eq. p7l39p with integrating Q!s(m^) over in the BS equation Doing so, we accounted for the 

analytical properties of all terms in the integrands of the BS equations, which had not been done in the preceding 
calculations. The further simplifications of aj-^^ depend on the integration region over /3. For example, if essentially 
/? ~ 1, Eq. ([M)) converts into the well-known expression a's-^^ « asik"^). The effective coupling ([57)1 explicitly depends 
on the value of n whereas a^-^^ in Eq. does not contain /i and looks /i- independent. However, it also depends on 
/i, though implicitly, through Eq. ([55]) . 

The expression for a^^-^ in Eq. pjp incorporates the contributions containing tt. They are originated by Eq. ([9]) where 
the analytical properties of as are respected. Those contributions are certainly beyond the logarithmic accuracy which 
is typical for the small-a; evolution equations, so there could be a conclusion made that such contributions should be 
dropped since very beginning, right in Eq. As the matter of fact, it cannot be done because the BS equations 
for the DIS structure functions and parton distributions contain contributions of the graph in Fig. 2 with the s-cut. 
Obviously, such a cut involves Sas(m^), i.e. in -terms in Eq. ((9]). So, neglecting the in -terms before the integration 
over automatically would mean accounting for the pole contribution psp only, i.e. would lead to fixing as at the 
fi -scale. However, Eqs. (|381 139)) show that after integrating over the contributions with tt can be neglected at 
large values of /i. 
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FIG. 1; The right-hand side of Eq. © 




FIG. 2: The integral contribution in Eq. ([TTJ 
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FIG. 3; Singularities of / given by Eqs. (IT4l- [T7Tl 
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FIG. 4: The integration contour C for calculating Ic 



